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Abstract 

The NJL model is formulated on the light cone. Using the 1/N 
expansion to solve the fermionic constraint, an effective 4-fermi la- 
grangian is derived. Pionic properties are investigated using this la- 
grangian, and the formal equivalence to the equal-time formulation 
is demonstrated. Two regularization schemes in terms of light-cone 
variables are discussed: An extension of the 'invariant mass cut-off 
scheme', and a transverse momentum cut-off scheme. It is shown that 
the first one is equivalent to the covariant 3-momentum cut-off (disper- 
sion cut-off) scheme in the equal-time formulation. As an application 
the structure function of the pion is studied in both regularization 
schemes. 
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1 Introduction 



As an effective theory of QCD in tlie low energy region, tlie Nambu - Jona- 
Lasinio (NJL) model is a powerful tool to investigate the properties of 
hadrons p|. It exhibits the spontaneous breaking of chiral symmetry and 
the appearance of Goldstone bosons in a simple and clear way. Due to the 
simplicity of the interaction, the relativistic two- and three-body equations 
can be solved in the ladder approximation, which allows a consistent unified 
description of mesons and baryons [Q, |[. Recently, extensions have been 
made to include also mesonic degrees of freedom into the description of three- 
quark bound states and methods to incorporate the effects of confinement 
phenomenologically have been proposed P]. Concerning the ground state 
properties of baryons, the NJL model has been quite successful . 

Recently, the model has been used also to describe the structure functions 
of the pion and the nucleon measured in deep inelastic scattering. It 
is well known that in the Bjorken limit the structure functions are expressed 
as Fourier transforms of correlation functions along the light-cone (LC) [p!0|] . 
Therefore, for the evaluation of the structure functions in the Bjorken limit 
the use of LC variables is most natural []TT|. The interpretation of the struc- 
ture functions as LC momentum distributions of the constituents (partons) 
in the ground state of the hadron, however, is strictly speaking possible only 
in a framework where the LC correlation functions can be considered as equal 



'time' correlations. In the LC quantization procedure |]1T], |I2|, |T3[ this is in- 
deed the case, since in this approach the fields evolve in the 'time' variable 
= {x^ + x^)/\^ instead of the usual time variable Therefore, from the 
physical standpoint one expects that the most natural way to describe the 
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structure functions in the Bjorken limit consists of not only using LC vari- 
ables, but also LC quantization. The LC quantization, which was proposed 
by Dirac [jl4[ a long time ago, has also many other interesting features, like 
the 'vacuum triviality' []T3| and associated with this the fact that Fock space 
expansions of physical states have an exact meaning, without referring to a 
particular approximation scheme [0. 

In spite of these physically appealing features, LC treatments also have — 
mostly technical — disadvantages. For example, they conceal the underlying 
rotational invariance and complicate the treatment of angular momentum. 



and bring in new problems of how to regularize divergent integrals 0] . More- 
over, the presence of zero modes destroys the 'vacuum triviality' in the naive 
sense, and the proper treatment of these zero modes is actually very com- 



plicated and technically involved |T^. Therefore, in actual calculations of 
structure functions in the NJL model, it might eventually be more conve- 
nient to use the usual equal time (ET) quantization and merely perform a 
transformation to LC variables at some convenient stage. Nevertheless, it 
is an important task to show that such a procedure is indeed equivalent to 
the LC quantization method with all its physical merits. To investigate this 
equivalence for the case of the NJL model is the purpose of this work. 

The demonstration of the equivalence of ET and LC quantization for 
cr-model type theories was the subject of previous investigations |T^. For 
4-fermi theories like the NJL model, it is well known that the LC quanti- 
zation poses additional problems due to the nonlinearity of the fermionic 



constraint [T^, [T9| . To circumvent this problem, in some previous works JlQ 



this fermionic constraint equation was replaced by the corresponding equa- 
tion for free fields. In this work we will show a straight forward method 
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based on the expansion up to the next-to- leading order to solve the 
fermionic constraint. (For methods using boson mapping techniques see ref. 

We will show that the resulting effective lagrangian has the same form 
as obtained by using the free field equations to solve the fermionic constraint, 
but the 4-fermi coupling constant is replaced by an effective one. For phys- 
ical quantities, the changes due to this modified coupling constant cancel 
against those due to modified propagators, such that the theory is formally 
equivalent to the ET theory. We will show that there exist cut-off schemes 
where this equivalence holds even for the regularized theories. Namely, a 
suitably extended version of the "invariant mass cut-off" introduced in ref. 
[ PT[] will be shown to be equivalent to the ET treatment based on the co- 
variant 3-momentum cut-off (dispersion cut-off) [^. We will also discuss a 
recently proposed alternative regularization scheme which has no coun- 
terpart in the conventional ET treatment. As an application, we will discuss 
the pion structure function, and compare our results with the empirically 
derived valence quark distributions in the pion [Q. 

The rest of this paper is organized as follows: In sect. 2 we review some 
well known facts of the NJL model in the ET formulation based on the co- 
variant 3-momentum cut-off scheme, and discuss the pion structure function 
in the ET formulation. In sect. 3 we derive the effective LC lagrangian in the 
framework of the 1/N expansion, discuss the new Feynman rules based on the 
LC quantization, and show the formal equivalence to the ET theory. Two 
regularization schemes suitable for LC theories will be discussed. In sect. 
4 we will show our numerical results for the valence quark LC momentum 
distribution in the pion, and in sect. 5 we summarize our work. 
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2 NJL model in the ET quantization and the 
pion structure function 



2.1 The model 

The lagrangian of the flavor SU(2) NJL model in its original form is given 
by i 

C = C^ + Ci, (2.1) 

where Cq = — m)ilj is the free quark lagrangian with m the current u, d 
quark mass, and Cj is given by 

£, = G((V^^)'-(^(75r)^)'). (2.2) 

The gap equation is given as usual by 

M = m- 2G{^ij) =m + UzG j^,t^D(S{q)^''). (2.3) 

M is the constituent quark mass, S{q) = - — ^ ^ . 7° is the Feynman prop- 
agator for the constituent quark [], and (V'V^) denotes the quark condensate in 
the leading order (A^) of the expansion. The solution of the BS equation 
for the qq t-matrix in the pionic channel is 

t(AO = (7%5T^>(A;)(7°75T^.). (2.4) 

^To correspond to the conventions used in most papers on LC field theories (see for 
example ref. p]), we define the Feynman propagator via the T-product of ijj and V'^, 
and normalize spinors according to {p, s')u{p, s) — 2EpSs's- We also restrict om'selves 
to the lowest order in the 1/N expansion (direct terms only), and therefore the effective 
coupling constant in the scalar and pseudoscalar channels is given by G (instead of jfG). 
Throughout this paper, < tpTtl; > denotes the leading order (N) of the vacuum expectation 
value of ^Tip. 
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Here the reduced t-matrix is given by 



where 



^2^^4'"^ ' 'o-^viy / /o'-'Vi (2-6) 

is the famihar 'bubble graph'. The pion mass m-Tr is obtained as the pole of 

(ED- 

The Bethe-Salpeter (BS) wave function ^hifl) in the covariant normaliza- 
tion is defined via the residue of the connected Green function at the pion 

I2 



pole kQ = ujk = \lk + ml, i.e; 

S{q')S{q' - k)t{k)S{q - k)S{q) ^ — f ^ as fco - c.,, (2.7) 

which leads to 

<^>U?)=^?(%)7W(g-A:)); 1 = - ^. (2.8) 

Here g is the pion-quark coupling constant. The assignment of the external 
momenta in eq. ( |2.7| ) is shown in fig. ^ 

The pion decay constant is defined as usual via the matrix element of the 
axial vector current between the pion state \ki > and the vacuum: 

< 0|^(0)|7'^75V^(0)|A;^ >= (2.9) 

This amounts to calculating the familiar one-quark loop contribution to the 
pion decay, using the BS wave function eq. ( pTSf) at the qqn vertex. The 
results is 

zP/. = 3gJ-0^^trn{^'ri^S{q-kh'j,S{q)) (2.10) 



UgMk^' 1^ ^ ^ (e = ml). 

' (27r)4 (^2 _ M2 + ie) ((g - kf - + ie) 



(2.11) 



2.2 Quark LC momentum distribution 

In the ET quantization method, the valence quark LC momentum distribu- 
tion qv{x) in the pion as a function of the Bjorken variable x [|10| 



q^{x) = / dz-e'''''-'' {ki\ij{0)-f+ij{z-)\ki) (2.12) 



ATI- 
can be obtained from the diagram fig. where the external operator is 7"'" 
and the loop integral is performed with respect to g+ and with g_ = 
xk^ fixed f\. Using S{q)'y^'y~^ S{q) = —dS{q)/dq+ and performing a partial 
integration in g+, we get 

= 3z//^^±||^tr^(%)7V%)7%55(g-A;h%5)(2.13) 
{k' = mll (2.14) 



^dli{k,x) 2_ 2 



^ dk 
where 



n(A:, x) = J ^^0^trn {l'l,S{qh'^,S{q - A:)) (2.15) 



i2n) 



q- 



with ^ ~ J~- "^^^ quantity (|2.15|) is defined such that 



/oo 
dxU{k,x) = U{k). (2.16) 
-00 

To evaluate ( p. 15 ), we need the Feynman propagator in terms of LC variables: 



— + le 



■^Our conventions for LC variables are summarized in Appendix A. 



where 



with 



p2 — + ie 2p_ yp+ — Cp + p+ — — ie^ 



(2.18) 



(2.19) 



The Feynman propagator ( |2.17|) contains a contact term in the LC 'time' 
variable x"*" (conjugate momentum p+), and it is convenient to separate this 
term: 



S{p) 



, 7+7° , 7+7° 



-7" + 



(2.20) 



- M2 + ie ' ' 2p_ ' 2p_ 

where is the on-shell four momentum with LC components (p+ = ep,p^, Pj_] 
Inserting this into eq. (|2.15| ) and (|2.16|) , one obtains 



U{k,x) = U{k,x) + R{x), 



U{k) = U{k) + R, 



(2.21) 



where Il{k,x) [Il{k)j is obtained by replacing S ^ S in eq. (|2.15| ) (eq. 
roD). For k \ = we obtain 

60(x)e(l -x) f <fq^ ql + 



U{k,x) 



x{l — x) 



(27r)'^ ql + M^- Px(l - x) - ie 



U(k) = / dxU{k,x) 
Jo 



(2.22) 
(2.23) 



while R is due to the contact term and given by 
R = I dx R{x) - 



{2-nY \q--k- g_ 



(2.24) 

where V denotes the principal value with respect to the g_ integration. Due 
to eq. ( |2.2CI| ), S can be replaced by 5* in this expression. Using ( |2.18|) and 
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performing the g+ integral we obtain Q 



oo \x — 1 X — 1 X X J J (27r)^ 

(2.25) 

It should be noted that in obtaining the expressions (|2.22| ) and ( p.25|) , no 
shifts of integration variables have been performed. 
The distribution function ( 2.14|) then becomes 

= -9 = % / .o^^..,.. , o .. -X^' (2-26) 



dk^ ^ J i2n)^ (ql + - mlx{l - x)y 

and since = —dli{k)/dk'^ {k'^ = m^) (see eq. ( p.8|) and (|2.21|) ), the sum 
rules 

/ dxqy{x) = l] I dxxqy{x) = - (2.27) 
Jo Jo 2 

follow immediately. The physical content of these relations is that we have 
one valence quark in the pion which carries half of the LC momentum k-. 
Finally, from the following identities 

f d^q 1 ^ r , f d\x 6(q) 1 f d^q^ 

7:rT4T-n o . . .2 = 7; dq 



(27r)'(g2-^2 + ,,)2 27 ^ 7 (^27ifql+fJ.^ 2 J {2T,f ql+ ^j? 

(2.28) 

we note that: (i) One could rewrite the expression (|2.26|) for qv{x) in terms of 
a 4-dimensional integral, which agrees with the form given in previous works 
[|]. Such a procedure could be used in order to avoid the regularization in 
terms of LC coordinates, but we do not follow it in this paper, (ii) There is 



^We define 



dx f dx 

lim / — = =pi7r 



, a; — a ± ze r^qo J_ji x — a ±ie 

where the second equahty foUows from usual contour integration, taking into account the 
contribution from the semicircle. 
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no contribution from the vertex correction shown in fig.^ This graph, which 
arises from the structure of the constituent quark itself, gives a contribution 
oc 6{x) when inserted into the graph of fig.@. Due to the regularization 
procedure discussed below, this gives no contribution. 

2.3 Regularization 

In the present ET formulation, we will regularize the gap equation (|2.3|) , 
the bubble graph ( p.6| ) and the expression for the pion decay constant ( p.llj ) 
using the covariant 3-momentum cut-off (A3) as explained in ref. (This 
scheme will be denoted by '3M'.) Basically, this consists in integrating over 
go using the theorem of residues, and then cutting the integral by |qf| < A3. 
Since this procedure is not Lorentz invariant, however, for quantities with 
two-body intermediate states like (p^ ) and ( p.ll|) it is necessary to specify 



the frame where the integral is cut off, and then to 'boost' the result to a 
general frame. If we cut off the integral in the frame k = 0, the 'boost' for 
Lorentz scalar quantities simply means to replace fcg — ^ k'^ in the resulting 
expressions. In this way we obtain 



n3M(fc) = 48/ /^l^ f; = f dfi 



q|<A3 (27r)3 fc2 _ 4^2 _^_ ,g J^^j, ^ fc2 _ ^2 + 

(2.30) 

where Eg = A/M^~+~g^, and the subscript 3M indicates that these quantities 
refer to the 3-momentum cut-off procedure. In the second equality of eq. 
( |2.30|) we introduced /i^ = 4(M^-|-q^) to express 11 (A;) as a dispersion integral 
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with p(/i2) = ^ f 1 - ^ j , and A2 = 4(A2 + M'). 

The following important relation between the condensate and the bubble 
graph at A; = 

{^t/j) = Mn(0), (2.32) 

which can be derived from the chiral Ward identity for the pion decay vertex 
(see Appendix B), is satisfied in the 3M regularization scheme, i.e; {4'ip)3M = 
Mn3A/(0). Therefore the gap equation (|2.29|) and the pion pole condition 
1 + 2GIl{ml) = can be written as 

m 

l + 2Gn3M(0) = — (2.33) 
n3M(m^)-n3M(0) = <nml) = -^, (2.34) 

where /(A;^) was introduced in ( |2.31| ). Eq. ( |2.34D shows explicitly the Gold- 
stone boson nature of the pion. The GOR relation is then obtained as follows: 



Mm m\ ,,9 9x ''^ \ 

m) 

1 dlim?) 

with C = — — ^ ^ ■ In the last step we used the relation = — 12M^ 
/(m^j om^ 

xI(ml)/(l+mlC), which follows from (OTI) and 4 = (k^H^^)) o , • 

For quantities expressed in terms of LC variables, the 'invariant mass 
cut-off' introduced by Lepage and Brodsky [^] is very convenient since it 
preserves all kinematic symmetries. (We will denote this cut-off scheme by 
'LB'.) It is applicable to Feynman diagrams containing intermediate states 
with n > 2 particles propagating with respect to the LC 'time' variable . 
For example, the bubble diagram n(A;, x) involves a qq intermediate state 
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with invariant mass squared 2A;_ (e^ + Ck-q) = (gi + M^) —x). (We set 
k± = 0.) Requiring this to be smaller than leads to the 'LB condition' 

LB: ql < A^x{l - x) - — > Xi < x < X2, (2.36) 



where xi^2 ~ 2 \/l ~ AM'^/A'^^ . This prescription regulates the ultra- 
violett as well as the infrared divergencies. The resulting regularized expres- 
sions IlLB{,k,x), UiBik) and qv,LB{x) are then given by eqs. ( p.22| ), ( p.23| ) 
and ( p.26|) with the condition ( p.36|) . 



A crucial question of course is whether this LB regularization procedure 
is consistent with the 3M procedure used above. For this, we note that if we 
perform a variable transformation x with 

^ = ^1 + ;|) , iEq = ^jql + ql + M2) (2.37) 

in the expression for nLB(A;), we exactly obtain H^Mik) of eq. ( p.30| ) with 
A^ = 4 (A3 + M^). That is, the LB cut-off A is the same as the dispersion 
cut-off A introduced in ( p.30| ). Similarly, if we evaluate the expression ( p.llj ) 
for in terms of LC coordinates, 

= QgM ( dx I ^ , -, (2.38) 

•'^ ^ J J {27f)^ql + M^ -mlx{l-x)' ^ ' 

impose the LB condition ( |2.36| ) and transform a; — > gs as above, we arrive 



at the expression eq. ( p.31[ ). We therefore see that the 3M and LB cut- 



off procedures are consistent in the sense that (i) they lead to the same 
expression for /tt, and (ii) the quark distribution regularized with the LB 
cut-off involves the quantity YiLB{k,x) (see eq. ( |2.26| )), the x-integral over 
which coincides with the bubble graph Yl'^Mik), i.e; IlLBik) = U^Mik). 
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At this point the question arises whether it is possible to consistently 
regularize all integrals in terms of the LC coordinates from the beginning, 
without making recourse to the regularization in terms of the usual coordi- 
nates. For this purpose, the LB prescription has to be extended, and this 
will be discussed in subsect. 3.5. 

3 NJL model on the light cone 

In the ET quantization method used in the previous section, the introduc- 
tion of LC variables is merely a change of variables. In this section we will 
investigate if we can obtain the same physical results in the framework of LC 
quantization. For this purpose, the main task is first to derive the form of the 
classical lagrangian in terms of the dynamical degrees of freedom, from which 
the hamiltonian and the canonical quantization conditions can be derived. 

It is well known that in LC field theories not all components of the fermion 
field tp are independent degrees of freedom |TT|. Namely, if we introduce 
projection operators 

A± = (3.1) 

and split the fermion field according to ip = A+t/^ + K^il) = + ip^, the 
component ip- is not a dynamical degree of freedom since there is no 'time' 
derivative 9+ acting on it in the lagrangian (see eq. (|3.4| ) below). Therefore, 
ip- should be eliminated by using the equation of motion before the theory is 
quantized. This can be easily done for theories which contain only fermionic 
bilinear terms, like free field theories or Yukawa theories 1113 . For 4-fermi the- 



ories, however, the 'fermionic constraint equation' cannot be solved exactly 
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In subsect. 3.1, we will use the 1/A^ expansion to derive the effective LC 
lagrangian for ip^ . Before plunging into the formalism, let us summarize here 
the result: The lagrangian is expressed in eq. (p.20|) , which shows that the 
interaction term has the form obtained by naively using the free field equa- 
tions to eliminate ip- (see eqs. (|3.9| )), but there appears an effective 4-fermi 
coupling constant G (see eq. ( |3.21| )) instead of the original G. In subsect. 
3.2 we discuss the Feynman rules which follow from this effective lagrangian, 
and in subsects. 3.3 and 3.4 we re-derive the qq t-matrix in the pionic chan- 
nel, the pion decay constant and the quark LC momentum distribution in 
the pion. There we will show that the results for these quantities are for- 
mally the same as obtained in the ET theory. In subsect. 3.5 we discuss two 
possible regularization procedures in terms of LC coordinates: An extended 
version of the LB scheme, and an alternative method which avoids the 
introduction of a cut-off for the longitudinal momentum components. 

3.1 Effective LC lagrangian in the 1/N expansion 

G 

To define the 1/A^ expansion, we write as usual G = — with Go independent 

of A^. Let us first discuss the 0{N) of the terms in the lagrangian ( p.l| ). 

Splitting ijjijj =< ipi/j > + where < ... > denotes a contraction with 

respect to the fermi fields, we have < ipip >oc and iptp -?j oc \/N. 

Therefore the lagrangian ( p.l| ) contains a c-number part of 0{N), a 1-body 

part of order 0{^/N), and a 2-body part of order 0{1). Our aim here is to 

give a similar decomposition of the lagrangian in the LC theory, i.e; to derive 

the 1-body part to 0{\^) and the 2-body part to 0(1). 

^This follows from < (: ip^jj :) >cx N. A similar consideration shows that in the 
decomposition (: tpip :) =< (: 1/;^/; :) > + (:(: tpip :) : J both terms are of 0{N). 
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To formulate the expansion in the LC theory, it is convenient to 
introduce bosonic auxihary fields a and tt. The lagrangian (|2.1j ) is equivalent 

to 

C = ^{^ - M)^ - ^ (a^ + TT^) , (3.2) 

where 

G 

M = m + 2^{a + i7r- T75) . (3.3) 

We then rewrite this lagrangian in terms of LC variables, making use of the 
projection operators ( p.l|) . Some useful relations for quark bilinears are given 
in Appendix A, and we obtain 

- i^l (ia^ ■ d± + -i^m) ip_ - i)l (ict^ ■ d± + 7°m) (3.4) 

Since we introduced the bosonic auxiliary fields, the fermionic constraint 
equation dC/dtpl^ = can be solved explicitly: f\ 

iP_ = -j}— (m^ ■ + -f'^M) (3.5) 



and inserting this into eq. (|3.4| ) we obtain 



[tct^ ■ + 7°m) ^ {ic^i. ■ + 7°m) (3.6) 



^We use the following definition of — ^: 



with 6(2;) = e(x) - Q{-x). 
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We now separate the leading order (N) of the 'c-number part' of the a 
field by writing a = a^N + a. Here a contains the leading order (V^) of 
the 'fluctuation part', as well as the next-to-leading order (A^'') of both the 
c-number and the fluctuation part. From the form of the lagrangian ( |3.6D it 
is clear that up to the order which concerns us here (that is, 0{\/N) for the 
1-body part, 0{N^) for the 2-body part) it is sufficient to consider a (and 
also tt) as the fluctuation piece of C(-\/iV). We define the constituent quark 
mass as 

M = m + 2Goao. (3.7) 
To split the lagrangian into 'free' and 'interaction' parts, it is convenient to 



follow ref. 1 12] and define the fermion field 

^ = ^+ + (3.8) 

with 

vi>+ = vl/_ = {ta^ . + 7°m) V+. (3.9) 

In terms of this field the lagrangian (|3.6|) can be written as 



C = Co-2^^{a + t7T-T^,)^-^{a' + n^)-2Goaoa 

- 2(^^) ^(a + «7r-r75)7+T^('^ + ^7r-r75)^ (3.10) 



with 



£o = ^(z^-M)^ 

V2i;iid+^+ - -^^t {la^ ■ d± + 7°m) {ict^ ■ d± + -f^'M) 

(3.11) 
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We note that up to the required order it is sufficient to replace the two fermi 
fields in the last term of (|3.10|) by a contraction. To exhibit the dependence 
of the lagrangian on the bosonic auxiliary fields more explicitly, we express 
it in terms of the combined field $ = {a, tt) (i.e; = a, = tt*) and 
the corresponding vertex operator F = (l,i75T) as 

C = Co- ^•I'^^^A^^'I'^^) - <l>(")5„, (3.12) 

where we introduced 

Aa, = 2^5.,+ (2^)'(^r.7+-i^r,vl/) (3.13) 

1 0— 

= 2(5,0^0^0 + 2^^1,^. (3.14) 

We now eliminate the auxiliary boson fields using dC/d^ = — yl$ — 5 = 0, 
which gives formally 

<i> = -A-^B; C = Co + -B'^A-^B. (3.15) 

The matrix operator A introduced above is defined by 

= I (a^, <I>(^) + <l>('^) A,,) . (3.16) 

Due to parity and/or isospin invariance, A^/s oc 6a/3, and we therefore define 

Aa,^6^,2^ (i + 2^f), (3.17) 



The quantity A involves = ^ acting to the right, while A involves acting to 



the left. The symmetrized form of A in (3. IE) appears naturally if we take the functional 
derivative of the action J d'^x£{x) with respect to using the definition of 1/d-. We 

note that also in the lagrangian eq. ( |3.12 ), Aaf3 can be replaced by Aaf3, since this leaves 
the action unchanged. 
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where the quantity 

F= i I <^7+^^^ > + <^7+^^$ > I (3.18) 
2 V id- i d- I 

is of O(A^). 

(To is determined by the requirement < $ >= —A~^ < B >= 0. Due to 
( p.l4| ) and (|3.7|), this requirement < B >= leads to the gap equation 



M = m-2G'<^^>. (3.19) 

Since in this leading order (A^) we have < \1'\1' >=< ipip >, this is formally 
the same gap equation as obtained in the ET theory, eq. ( p.3|) . 

Since < B >= 0, the interaction term in the lagrangian ( p.l5| ) becomes 
B :) A-^ {■ B :). We therefore arrive at the final form of the LC la- 
grangian: w 

£ = £0 + G (^(: :)' - (: ^75T^ :)') , (3.20) 



where Cq of eq. ( |3.11| ) contains the 1-body part of 0{yN)^ and the 2-body 



part is of order A^" and involves the effective coupling constant 

with G = Gq/N. One should note that the lagrangian (|3.2CI|) contains only 
ip^ as a dynamical component, since we defined the field \E' in terms of 
in eqs. ( ^.81 ), ( |3.9| ). Our result is that the effective interaction lagrangian 
in the LC theory up to 0{N^) can be obtained from eq. \2. i| j by using the 
free field equations to eliminate ip^ , but the coupling constant G is replaced 



^Our notation used in cq. (3.2C) anticipates the fact derived in the next subsection 



that in coordinate space the quantity G (x 5(x ~ y ), i.e; that the interaction lagrangian 



is local although F of eq. (3- IS) is apparently non-local. 
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by the effective coupling constant G of eq. ( ^.2]\ ). (The explicit expression 



for the quantity F will be given in eq. ( p. 34 ) below. 



Once we have derived the classical lagrangian for the quantization 
proceeds as usual by imposing anti-commutation relations between ip^ and 
n+ = v^^^+ at equal LC times x~^. The mode expansion of the free field 
reads f\ 

with d^p = dp-d'^p± , p = (p_,pj_), p ■ X = p+x"*" + p^x~ — ■ x±, and 
was defined below ( p.20|) . Using the mode expansion given above, one 



^+(x) = J f^^^^ L_ Y: [bp,suAp, s)e-^^-- + 4,,t;+(p, s)e^^-] (3.22) 



obtains for the free fermion propagator 



S++(p) = F.T.(-2<0|t(^+4)|0>) 



(3.23) 



where T orders the operators with respect to x^ . The quantity l/{p^ — + 
ie) has been given in terms of LC variables in eq. (|2.18| ). 



3.2 Feynman rules, gap equation and effective cou- 
pling constant 

To calculate a fermionic Green function 

<0\f{ij+{xi)...tPi{xn))\0>, (3.24) 

one could formulate the Feynman rules directly in terms of the field ?/'+, 
using the lagrangian ( |3.2CI| ) and the definition of \1/ in terms of given in 



^Properties of the LC spinors u± — A±u and v± = A±v are summarized in Appendix 

A. 
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eq. ( |3.8| ). In this method, one uses the propagator eq. ( p.23| ) for internal 
fermion hnes, and 4-fermi vertices of the following form, referring e.g; to the 
pionic qq channel, 

-2{GXl{q') Xl{q) 
^ -2iG[^,^'n{Z{q')-Z{q' -k))] [^,^\,{Z{q - k) - Z{q))\ , 

(3.25) 

where the external momenta are assigned as in fig. |l], and 

Z{p) = (oc^ ■ + 7OM) (3.26) 

A much more convenient method, however, is first to calculate the Green 
function 

<0|f (^(xi)...^t(a;^)^|0>, (3.27) 

and then to project out ( ^.24| ) by applying the projection operator A+ to each 
external fermion leg. In this method, which will be used in the next section 
to obtain the qq t-matrix in the pionic channel, one uses the propagator 

S(p) = F.T. (-t < OIT f |0 >) = , ^t.!^ 7° = S(p) - 

(3.28) 

for each internal fermion line, and (constant) 4-fermi vertices which differ 
from those in the ET treatment of the NJL model only by the replacement 
G ^ G. The form of S{p) given in eq. (|3.28D can be easily derived as follows 



1^ : Splitting S{p) into the four components = +, — ), where 5*++ 



is given by eq. ( |3.23| ), we have, using the form of "if as given by eq. (|3.8[) and 
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the Feynman propagator ( |2.17] ) (see also Appendix A): 



5+_(p) = 5++(p)Z(p) = A+S(p)A_ 
S^4p) = Z(p)5++(p) = A.5(p)A+ 

S^4p) = Z(p)S++(j9)Z(p) = A_5(p)A_ - . (3.29) 

V2 p- 



/^From these relations one easily obtains eq.( |3.28| ). This propagator S{p) has 
been introduced already in eq. (|2.20|) by separating the contact term from 
the Feynman propagator S{p). 

Clearly, if the ET and the LC treatments of the NJL model are equivalent, 
the contributions due to additional term —j'^'j^ /2p_ in the propagator and 
the modifications due to the replacement G ^ G in the vertices should 
'cancel' for physical observables. These are the same kinds of cancellations 
as observed for the a-model in ref. ||12|| . 



Using the propagator (|3.28| ), it is now easy to derive the explicit forms of 
the gap equation and the quantity F, which determines the effective coupling 
constant (see eqs. ( p. 18 ) and ( |3.21 )). For the gap equation we obtain from 



(13) and ( CT) 



M = m- 2G{^^) =m + l2iG j^^it^D {Siqh'') , (3.30) 

and since from eq. ( |3.28|) we have trD (^S{q)'y^^ = trn {S{q)'y^), we see that 
the gap equation is formally the same as in the ET theory. Using ( |2.18| ) and 
performing the g+ integral we obtain 

M-m = 48.GM//^^— -L = 24GMr^/^. (3.31) 

' (27r)4 g2 _ M2 + ie Jo g_ J {2nf ^ ^ 

Although it seems from this expression that in the case m = only the trivial 
solution (M = 0) exists, we have to note that the theory has not yet been 
regularized. (See subsect. 3.5.) 
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Let us now turn to the quantity F of eg. ( P-18| ). Using the definition of 
given earher, we have for the quantity A^p = Sa/sa of eq. (|3.17|) 



a(z-) = 2G 



6{z-) + 3Ge{z-) tro (S(^-)7V - S(-;2-)7V)] (3.32) 



where z = x — y ■ The Fourier transform of this is 



a(fc_) = 2G [l-QiCj -7r\,trn (S(g)7% 



(27r)* V ^-^^ ' ' J\q_ + k_ q_-k^JJ 

(3.33) 

Again we note that in this expression S{p) can be replaced by tlie Feynman 
propagator S{p) without any modification. Performing the g+ integral using 
( |2.18| ) and introducing g_ = k^x we see that a (or F) is actually a constant, 



and so is G of eg. ( p. 211) . The expression for F becomes 



oo 



3.3 Pionic t-matrix 

According to the discussion at the beginning of the previous subsection, to ob- 
tain the qq t-matrix in the LC formulation, we first determine the t-matrix as- 



>. 



sociated with the amputated connected Green function < 0|T (^\1/\I'\E'^\E'^ j |0 
In the ladder approximation this is obtained by summing the bubble graphs 
Il{k), which have been introduced already in sect. 2, in the same way as in 
the ET theory as 

t(fc) = (7°75r.)r(A:)(7%5r.) (3.35) 



with 



-2iG 

rik) = , , 3.36) 

n(A;) = 6z|^tr^(7°75%)7%5%-fc)). (3.37) 
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The t-matrix associated with the amputated connected Green function 

< 0\T |0 > is then simply obtained by projection (see fig. |l| 

for the assignment of the external momenta): 

S++{q')S++{q' -k) i S++{q - k)S++{q) 
= A+ (S{q'h'^,nS{q' - k)) A+ T{k) (S{q - kh'^,nS{q)) A+ 



(3.38) 

The result of this projection is t = Xl{q') r{k) Xl{q), i.e; the vertices (7°75Tj) 
are replaced by the momentum dependent vertices XI of eq.( p.25| ). Of course, 
the pole structure is not affected by this projection. 

Using eq.( p.21| ), we can rewrite the reduced t-matrix ( p.36| ) as 

-2iG 



(3.39) 



1 + 2G (n(fc) + f) 

Comparing with eq. ( p.5| ), we see that we should have the relation 

I{{k) = Ii{k) + F (3.40) 

with F given by ( |3.34] ), in order that the reduced t-matrices, and in particular 
the pole parts, in ET and LC theories are the same. On the other hand, by 
explicit evaluation we found that n(A;) = n(A;) + R (see eq.( p.21|) ), where R 
is given by (|2.25| ). Formally, the relation 

F = R (3.41) 

holds, but to establish it from ( p. 251) and (|3.34|) one has to introduce variable 
shifts in divergent integrals. 

The Goldstone nature of the pion in the present LC treatment follows by 
noting that the relation ( p.32|) , which now reads 



< >= M (n(0) + R 
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(3.42) 



is formally satisfied. The explicit verification by using ( |3.31| ), ( |2.22| ) and 



( 2.25 ), however, again involves variable shifts in divergent integrals. If one 



nevertheless assumes the validity of ( p.42|) and (p.41|) , one can write the gap 



equation ( ^.3(J ) and the pion pole condition 1 + 2G n(m^) = as follows: 



l + 2Gn(0) = - (3.43) 
n(m^)-n(0) = = -1^^ (3.44) 



where /(A;^) was introduced in ( p.31| ). Here we used ( p.21| ) and (p.32|) to 
express the gap equation in terms of G and 11(0), and introduced 

_ m G 

m = = m-. (3.45) 



Eq. ( |3.44 ) shows the Goldstone nature of the pion. We note the the GOR 



relation is satisfied with the current quark mass m (not m). Its derivation 
follows the same steps as outlined in eq. ( p.35| ). 



3.4 Bound state matrix elements and pionic LC wave 
function 

In contrast to the BS wave function introduced in sect. 2, a 'LC wave function' 
usually refers to the solution of a Schrodinger-type (or, more precisely, TDA- 
type) equation which can be derived directly in the hamiltonian formalism 
(see Appendix C). Here we want to show how this kind of wave function can 
also be obtained by using the projection methods introduced in the previous 
sections. 

In the same way as in the ET theory, we first define a BS wave function 
$ in the LC theory by 

— — — — i W,(a'W,Ua) 

S{q')S{q' - k) t S{q - k)S{q) ^ -—-Aiil-^ as k+ ^ e,. (3.46) 
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This wave function ^k{(l) is given by ( p.8| ) with S S. The BS wave function 
for the dynamical (+) fermion components ($) is defined by (see eq. (|3.38|) ) 

(3.47) 

It is clear that 



^l{q) = A+$fc(g)A+ = g[^\^S{qh'^,nS{q - k)A^j (3.48) 
Due to = 7~''7°A+ = and (p.28|), the propagator S can be replaced 



by the Feynman propagator S in this expression, and $ agrees with the ++ 



component of the BS wave function ([278 



The calculation of bound state matrix elements in the LC theory can 
then proceed by using these BS wave functions in the same way as in the 
ET theory. As in the case of the Green functions, one has two choices: One 
is to use the wave function $ together with the vertices and propagators 
defined in terms of the field and the other, more convenient method is to 
use the wave function $ together with the vertices and propagators defined 
for the field In Appendix D it is shown that the calculation of the pion 
decay constant and the valence quark distribution in the pion along these 
lines leads to the same results as in the ET theory, namely eq. (|2.11| ) and eq. 



Another method to calculate matrix elements involving bound state like 
those in ( p.9|) and ( |2.12|) , which is more specific to LC theories, is to use the 



Fock space expansion of the pion state \ki >. Let us explain this approach 
here in order to make contact to the concept of Schrodinger-type LC wave 



functions. It is well known |25| that, for interactions which do not depend on 
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the 'time' component of the relative momentum (g+ here), the Schrodinger- 
type (equal 'time') wave function is obtained by integrating the BS wave 
function over g+ : 

tf )(,) = / p±m,) 

J 2m 

= ^ / ^ (A+:S(g)7Vr.:S(g - A;)A+) . (3.49) 

As shown in Appendix C, this 'LC wave function' satisfies the same Schrodinger- 
type equation as derived directly in the hamiltonian formalism. If we use 
the propagator ( ^.28[ ) with the pole part given in eq. ( |2.18| ) in this expres- 
sion and consider the case > 0, we get the restrictions g_ > and 
— < 0, and therefore the LC wave function (|3.49| ) is non-zero only 



for < X < 1, where x = q-/k-. To go from the Dirac to the spin repre- 



sentation, we multiply eq. (|3.49| ) from left by [-\/2g_] ^ u+{<l, Si) and from 



— 1 /2 

right by (^V^{k^ — g_)^ Vj^{k — q, 82)- Multiplying also the corresponding 
isospin (t) and color (c) wave functions, and denoting a = {s,t,c), we obtain 
after a simple calculation 



$^(5) (a) = -i- ^ I (y\tyA5 



cl,c2 



X MUq,Si)75(7±-q± + M)t;+(fc-q,S2). (3.50) 



The conventional normalization 
72 



l^.fd^Y. \^&\aM- = k~x,q^)\' = 1 (3.51) 

1 

can be checked from eq. ( ^.5U| ), using 1/g = —dU/dk"^ ^ (/c^ = m^) with 
n(A;) given by ( p.23|) . The Fock space expansion of the pionic bound state in 
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the TDA is then given by 



H >= V2 f~ dq. [ ^Z\aM%,aA-nJ^ > ■ (3-52) 

"'0 aia2 



If we note that V(0)7+V('2") = y2^+(0)V'+ 

{z ) in eg. ( p.l2|) , we can use the 



mode expansion (|3.22| ) (supplemented by the isospin and color parts) and 
the Fock expansion (|3.52|) to obtain after a simple calculation 



= / ^ E l^ZUl- = k-X,Q±)\'- (3.53) 

The evaluation of the sum then leads to the same result as obtained in the 
ET theory, namely eq. (|2.26| ). 



In a similar way, by noting that the + component of the axial vector 

current is given by —^il'+l5''"^+y we obtain from eq. ( p^.9D the expression 
V 2 

(|2.38|) (yet unregularized) . 

3.5 Regularization schemes in LC theory 

In the previous subsections we have shown that the NJL model on the LC 
gives formally the same results for the pion decay constant, the gap equation, 
the pion pole condition and the quark distribution in the pion as the ET 
theory. In sect. 2 we noted that the bubble graphs in the LC theory (like 
Il{k, x) and n(/c)) can be regularized in a covariant way consistently with the 
3M regularization scheme in the ET theory. To give consistent regularization 
prescriptions for tadpole graphs in the LC theory, like the quark condensate 
or the quantity F of eq. ( p.34| ), turns out to be a very difficult task which we 
do not attempt to solve explicitly in this work. We will instead make use of 
the general relation (|3.42| ). We have shown that this relation holds formally 



in the LC theory, if one performs variable shifts in divergent longitudinal 
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momentum (g_) integrals. Here we will assume that ( p.42| ) and also ( p.41j ) 
are valid even in the regularized theory. This method of defining the quark 
condensate in LC theories in terms of correlation functions has been used in 
several recent works [ |26| | . 

Following this method, we will give two possible regularization schemes 
below. The first one ('extended LB scheme') is equivalent to the 3M cut-off 
scheme in the ET formulation. The second one ('transverse cut-off scheme'), 



which was proposed in ref. |]23|, avoids the introduction of a cut-off for the 



longitudinal momentum component. 
3.5.1 Extended LB scheme 



In addition to the LB condition ( 2.36 ) for the 2-point functions, in the 'ex- 



tended LB scheme' one assumes the validity of eqs. ( |3.42| ) and ( |3.41| ). Then 
the gap equation and the pion pole condition assume the forms ( |3.43|) and 
(|3.44D with n Hlb- In particular, the gap equation takes the explicit form 



dx f d'^q± m 
LB X J LB (27r)^ ~ M'' 



1-24G i - LP^ = ^, (3.54) 



where the LB condition is given by ( |2.36| ). Performing the variable transfor- 



mation (|2.37|) , we see that this is just the gap equation in the 3M regulariza- 
tion scheme (eq. ( |2.29| )) with G ^ G and m — m. 

Let us then compare a treatment based on the 3M regularization in terms 
of the usual variables with the 'extended LB regularization' in terms of LC 
variables. We treat the constituent quark mass M as a free parameter, choos- 
ing the same value in both regularization schemes. Since the expressions for 
/tt are the same in the two treatments, we obtain the same value for A in 
both schemes by requiring that the experimental value for is reproduced. 
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Then the 4-fermi couphng constant is determined by requiring that the ex- 
perimental pion mass is reproduced. In the 3M scheme this is expressed as 
1 + 2GIlsMi'rnl) = 0, while in the LB scheme we have 1 + 2G TlLB{ml) = 0. 
Since Ulb = ^sm, this gives G^m = Glb- Then the gap equation can be used 
in order to determine the value of the current quark mass. As is clear from 
the above discussion following eq. ( p.54| ), this gives fniB = ^sm- Therefore, 
the current quark masses (m) are different in the two schemes, but also the 
quark condensates are different: /^From ( |3.42| ) and ( |3.41| ), as well as from 
= we have 

{^^P)lb = imsM + MFLB (3.55) 
rriLB = (3.56) 

As a consequence the GOR relation holds in both schemes: In the 3M scheme 



we have eq. ( 2.35 ) with m = m^M, and following the same steps in the 



extended LB scheme we obtain 

- m,Bm,s = flml (1 - ^) (1 + mlG) . (3.57) 



M J 

This 'invariance' of the product of the current quark mass and the quark 
condensate holds also for m{H\ipip\H) , where \H) denotes a hadronic state, 
and means that in actual calculations the regularized form Flb and the 

value of m in the LB scheme are never needed: From the Feynman-Hellman 

- dMn dM 

theorem we have {H\ih'ih\H) = — — . The second factor in the latter 

^ ' ' ' dM dm 

expression represents the 'vacuum contribution', and its product with m can 
be written in the form (see Appendix B) 

dM m m 

9^ " 1 - 2Gn,(0) " l-2Gn<,(0)' 
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where Ho- (A;) is the bubble graph in the a meson channel, and T{fj{k) = 
Ilcr{k) — F. /^From our above discussion it is then clear that the product 
( |3.58| ), and therefore also the V-term' m{H\ipilj\H) , takes the same value in 
the ET treatment based on the 3M cut-off and in the LC theory based on 
the extended LB cut-off. 

The above discussion shows that, although the parameters G and m in the 
LC theory based on the extended LB cut-off scheme are different from those 
of the ET theory based on the 3M cut-off scheme, the physical predictions 
are the same. In this sense, the parameters G and fn introduced in the 
LC theory by eqs. (|3.2f|) and ( p.45|) can be considered as 'renormalized' 



parameters, which take the same values as G and m in the ET theory. (The 
actual value of the 'renormalization constant' 1 + 2GF is never needed.) 

We conclude that the extended LB scheme is physically equivalent to the 
3M scheme. 

3.5.2 Transverse (TR) cut-off scheme 



In this scheme , one notes that the divergence of the bubble graph 11 with 
respect to the longitudinal momentum integration is logarithmic (see ( p.22| ) 



and ( 2.23 )), and therefore one subtraction is sufficient to make it finite. As 
a consequence, the following equations are already finite with respect to the 
longitudinal momentum integration: (i) The pion pole condition written in 



the form ( p.44| ), (ii) the pionic t-matrix ( |3.36| ) written in the form 



nk) = (3.59) 

n(fc) — n(m^) 

(iii) the pion decay constant, and (iv) the quark distribution function. The 
TR scheme now consists in introducing a transverse momentum cut-off Atr 
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in these equations, while the remaining divergent equations (the gap equation 
(|3.43|) , and eq. ( p.34| ) determining the quantity F) are not used exphcitly. 

Note that in this method the gap equation is used only to formally rewrite 
the pion pole condition into the form ( |3.44 ), which is finite with respect to 
the longitudinal momentum integration. It is easy to show that this method 
works also for the sigma meson channel, i.e; the gap equation can be used 
to rewrite the equation for the sigma meson mass into a form which is finite 
with respect to the longitudinal momentum integration. It is a disadvantage 
of this method, however, that it works strictly only within the present 
expansion scheme: The gap equation and the eigenvalue equation for the 
meson mass must be of the same order in For the lagrangian ( |2.1| ), 

this is true only for the pion and sigma meson, since the interaction in other 
mesonic channels is due to the exchange term and suppressed by one power of 
It is, nevertheless, possible to apply the TR scheme also to these cases 
by performing a kind of mass renormalization, that is to use the eigenvalue 
equation for the meson mass only in a formal sense to rewrite the t-matrix 
into a renormalized form analogous to (p.59|) , which is finite with respect to 
the longitudinal momentum integration. Such a renormalization procedure, 
in which the gap equation and the eigenvalue equations for the hadron masses 
are used only in a formal sense, is quite different in spirit from the common 
usage of the NJL model. 
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4 Numerical results for the quark LC mo- 
mentum distribution in the pion 



In this section we present numerical results for the quark LC momentum dis- 
tribution ( |2.26| ) regularized by using the two regularization schemes discussed 
in subsect. 3.5, i.e; the extended LB scheme and the TR scheme. We choose 
the constituent quark mass M = SOOMeV in both cases. Then by requiring 
that ( p.31|) (or (|2.38| ) with the condition (|2.36| )) reproduces the experimental 
value = 93MeV, we obtain for the LB cut-off A = lAlGeV (correspond- 
ing to a 3-momentum cut-off of A3 = {).Q7GeV). ^ Similarly, for the TR 
scheme we use ( p.38|) with |qr_|_| < Ktr to obtain Ktr = 0.486*6^. The value 
of the pion-quark coupling constant, defined at k"^ = = (140MeV^)^, be- 
comes g = 3.15 (3.16) in the LB (TR) scheme. (These values are close to each 
other, since in the chiral limit we have g = M//^ = 3.23 in both schemes.) 

The resulting distributions multiplied by x are shown by the dashed lines 
in fig. ^ for the LB scheme, and in fig. |^ for the TR scheme. Due to the 
condition ( |2.36D in the LB regularization, the distribution qy{x) is nonzero 
only in the interval Xi < x < X2 and tends to zero smoothly at the end 
points Xi^2- In the TR regularization, on the contrary, g„(x) extends over 
the whole interval < x < 1 and is nonzero at the end points. In this case, 
due to the smallness of the pion mass, qv{x) shows only little variation as a 
function of x, which leads to the almost linear behaviour of xqy{x) shown in 
fig. 1^. This is very similar to the result of a previous calculation using 
the Pauli-Villars regularization based on Feynman integral representation of 
qv{x) (see the discussion below eq.( p.28|) ). 

^For the LB scheme, we obtain Gsm — Glb = i-bdCeV^^ from the pion pole condition 
(using ttxtt = 140MeV), and m^M — ttilb = 5.3MeV from the gap equation. 
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To connect these results with experiment, we assume as usual |T0[ that 
our NJL results are valid at some low-energy scale and perform the 
evolution using the DGLAP equations ^\ up to the value = iGeV"^, 



where parametrized fits to the experimental results are available EM. We 



used the computer code of ref. |^ to solve the DGLAP equations up to the 
next-to-leading order, using the parameters Aqcd = 0.25GeV and Nf = 3. 
The value of Qq is considered here as a parameter to be adjusted such as to 
reasonably reproduce the experimental distribution. In both regularization 
schemes the value obtained in this way is Ql = O.lSGeV^^, i.e; Qq is of the 
order of the constituent quark mass. 

The results of the evolution are shown by the solid lines in figs. ^ and 
^, and are compared to the empirical distribution given in ref. p4| (dotted 
lines). The fact that in the LB scheme the theoretical distribution is too low 
for high X is connected to the LB condition (|2.36|) which gives a vanishing 
input distribution for high x (see the dashed line in fig. 4). The use of smaller 
values for Ql would lead to better agreement around the peak position, but 
to a larger discrepancy for higher x. A similar situation is encountered also in 
the case of the nucleon structure function . On the other hand, the result 
based on the TR scheme reproduces the empirical distribution almost exactly. 
Therefore, from the phenomenological point of view, a regularization scheme 
which does not restrict the longitudinal momenta seems to be favorable. 



Since the computer code used for the evolution |28 requires an input distribution 
xqv{x) which vanishes at x = 1, we artificially modified it for x very close to 1 such that 
it goes like (1 — x)" with some power n. {n = 10 was used in the actual calculation.) 
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5 Summary and conclusions 



In this work wc investigated the NJL model in the LC quantization and its 
relation to the usual ET formalism, putting special emphasis on the applica- 
tion to the pion structure function. We first reviewed the NJL model and the 
calculation of the pion structure function in the ET formalism. We showed 
that the Lepage-Brodsky (LB) cut-off scheme for the quark LC momentum 
distribution can be applied consistently in connection with the covariant 3- 
momentum cut-off (3M) scheme for the gap equation and the qq polarization 
propagator. 

The introduction of LC variables in the ET formalism is simply a trans- 
formation of variables. As we discussed in the introduction, however, the 
interpretation of the structure function as LC momentum distributions of 
the constituents in the ground state target is most natural in the LC quan- 
tization procedure. Due to this fact, and other physically appealing features 
of the LC formalism, we formulated the LC quantization of the NJL model. 
The essential part was the solution of the fermionic constraint equation using 
the expansion up to the next-to-leading order. We showed that, up to 
this order, the interaction lagrangian on the LC for the dynamical compo- 
nent of the fermion field has the same form as obtained by using the free field 
equations to solve the fermionic constraint, but the original 4-fermi coupling 
constant is replaced by an effective (renormalized) one. The Feynman rules 
for the LC lagrangian are very similar to those for the ET lagrangian: The 
propagators for internal lines differ just by a contact term, and the 4-fermi 
vertices differ by a multiplying factor. For physical quantities these differ- 
ences formally cancel with each other, leading to the same results as in the ET 
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theory. We have demonstrated this equivalence for the eigenvalue equation 
determining the pion mass, the pion decay constant, and the pion structure 
function in the Bjorken limit. In the LC quantization, these quantities can 
be evaluated either by using the language of Feynman diagrams as described 
above, or by using the Fock space expansions based on the hamiltonian for- 
malism. Both methods lead to the same results. 

In this work, we avoided the difficulty of regularizing tadpole-like dia- 
grams in the LC theory by relating them to 2-point Green functions at zero 
momentum. We then formulated two possible regularization schemes: The 
first one was an extension of the "invariant mass cut-off scheme" ||21|. We 
have shown that this scheme leads to the same physical predictions as the 
covariant 3- momentum cut-off (dispersion cut-off) scheme in the ET formu- 
lation, provided that in both treatments the parameters are fixed such as to 
have the same values for M, /^r, and m-j^. The application of this scheme to 
problems others than studied in this paper, like the nucleon structure func- 
tion, is straight forward |2^. The second one ("transverse cut-off scheme") 
employs a cut-off only for the transverse momentum components. In this 
scheme, the gap equation is used only to formally rewrite the eigenvalue equa- 
tion for the pion and sigma mass into a form which is finite with respect to 
the longitudinal momentum integrations. We noted that this method works 
only within the present expansion scheme, and to apply it to other me- 
son channels or to the nucleon, one has in addition to perform a kind of mass 
renormalization. Such a procedure, however, departs significantly from the 
common usage of the NJL model. 

As an application, we evaluated the quark LC momentum distribution 
in the pion numerically and performed the evolution to compare the re- 
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suits with parametrizations based on experimental data. Since in the LB 
regularization scheme the input distribution is zero in an interval near the 
end points, the slope of the evolved distribution on the large x side is some- 
what too large compared with the empirical distribution. Nevertheless, a fair 
agreement can be achieved by assigning a low energy scale Qg — O.lSGel^^ to 
the NJL model, which is similar to the values used in previous investigations 
of the pion structure function. On the other hand, in the TR regularization 
scheme one can reproduce the empirical distribution almost exactly using 
the same value of Qo- This indicates that for phenomenological applications 
a regularization scheme which does not restrict the longitudinal momentum 
components seems to be preferable. 
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Appendices 



A LC variables and Dirac matrices 

In this Appendix we summarize our conventions for LC variables and Dirac 
matrices, and give some useful identities which have been used in the main 
text. We use the LC variables = (^^^ i ^^): ~ ^ ("^o -'- "^s)' 
x\_ = —x±i {i = 1, 2). Then we have for the Lorentz scalar product p ■ x = 
+ p^x^ — ■ x±. In the LC quantization procedure, x~^ is the 'time' 
variable and the 'space' variables. Similarly, p+ is the 'energy' 

variable and p- , pj^ are the 'momentum' variables. The dispersion relation 
for a free particle with mass M is p+ = Cp = + M^) /2p-, see eq. ( |2.19| ). 
For the derivatives we use the conventions da = d/dx"' (a = +, — , 1, 2). 

Analogous definitions hold for the LC Dirac matrices 7^,7^ terms of 
the 7^. It follows that ^ = 7"'' (9+ + 7"(9„ — 7_l ■ d±. Note that 7+7+ = 
7~7~ = 0. The projection operators A± have been introduced in eq. ( |3.1|) . 
Equivalent useful forms are A-t = :^7°7^ = ~ '^(^ ±7^7^). We 

have the following useful identities: 



7 A, 



A+7' 



,+ 



7+A 



A_7 







(A.l) 



•+ 



A±7" = 7' 




(A.2) 



For the quark bilinears it follows that 



V24^± 



(A.3) 




(A.4) 



^7± ■ dxip 



(A.5) 
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We also give some identities which were used in deriving eqs. ( p. 291 ) and 
(Pi: 



A±(^+M)A± = Ai(-7^.p^ + M) (A.6) 

A±(^+M)7°A± = y2A±p^ (A.7) 

A±(y+M)A^ = I^Pt (A.8) 

A±(^+M)7°A^ = ^7^(7±-p± + M) . (A.9) 

The Dirac spinors are normahzed according to ^^(p, s')u{p, s) = 2Ep6s's 
and v'^Ip, s')v{p, s) = 2Ep6s's- The spinors = A+m and f+ = A^v then 
satisfy the orthogonahty and completeness relations 

(p, s')u+ {p,s) = vl (p, s')v+ {p,s) = V2p^ 6ss' (A. 10) 

Y,u+{p,s)ul{p,s) = s)t;^.(p,s) = v/2A+p_ . (A.ll) 

s s 

B Relations between the condensate and 2- 
point functions 

We first derive relation ( p.32|) in the ET formalism from the chiral Ward 
identity for the pion decay vertex. If we define J^^ = ip'^^'^^Tiip and vr* = 
ij^Ti , the chiral Ward identity follows as usual from the PCAC relation 
df,f/ = 2m7r': 

d,{0\T{r/{x)7r^{x'))\0) 
= 2m(0|T (tt'{x)tt^{x')) |0) - 2i6ij6'-^\x' - x){0\^'^\0). (B.l) 



To lowest order in the correlation function (0|T {7t^{x)tt^{x')) |0) is given 
by the series of bubble graphs. Considering the Fourier transform of ( [B.l| ) 
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for k = 0, we therefore obtain in the leading order of 

= 2^m ^^^^L _2i(0|V^V|0). (B.2) 

l + 2Gn(0) ^ ' ^ ' 

Since the constituent quark mass is defined hj M = m — 2G{0\'ilJip\0), this 
relation can be written in the form (0|-?/'^/'|0) = Mn(0), which is eq. ( p.32|) . 
Next we derive eq. ( p.58D from the general relation 

d{ipip) 



dM 



-n,(0), (B.3) 



where T\.fj{k) is the bubble graph in the a meson channel. Eq. (|B.3|) can 
be derived most easily by using the functional integral. It is easy to check 
that this relation holds in the 3M regularization scheme to leading order in 
i.e; g^^*^ ~ — no-,3j\/(0). Then, in our extended LB scheme, differ- 
entiating (|3.55| ) with respect to M and requiring that the result is equal to 
— no-,LB(0) = — (no-,LB(0) — Flb) leads to the requirement that Fib is inde- 
pendent of M. By differentiating the gap equation ( p73D with respect to M 
we obtain 

dM m m m ,s 

dm l + 2G^{^P^|J) 1 - 2Gn,(0) l-2Gn,(0)' ^ ' ^ 

G m 1 

where we used = = — = = — . This proves eq. (B.58|). 

G m 1-2GF ^ ^ — ^ 



C Bound state equation 

In this Appendix we derive the bound state equation in the TDA for the LC 
wave function in the pionic qq channel, using (i) the TDA equation in the 
hamiltonian formalism, and (ii) the BS equation. 
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C.l Hamiltonian formalism 

The LC hamiltonian derived from the lagrangian eq. ( p.20| ) is given in terms 
of by 



V2id. 

_ r , „ f -1 

+ G 

+ cr channel, (C 
where we show only the interaction in the pionic channel explicitly. We 



V2td. \ ^ V2i d- 



now insert ( |3.22| ) into the above hamiltonian and apply it to the pion state 
( |3.52|) , keeping only those terms which are again of the form ( p.52|) . (This 



corresponds to the TDA, and means that in the interaction hamiltonian only 
terms of the form b'' S bd are kept.) The result is expressed as an eigenvalue 
equation for the 'energy' = (fc^ + m\)/2k_, i.e; H\k) = ek\k). The 
procedure is tedious but straight forward, and therefore we will give only the 
final result. One ends up with a TDA equation, where the kernel has both 
direct and exchange terms. The kernel for the direct term leads to a color 
factor 3, and in order to be consistent with the leading order in the 
expansion used throughout this paper we keep only the direct term. The 
result is 



''~^k~;: 2k4l-x) ) ^'^'--^^^ 



■ k,ai 

Jo , /o^n _ -rwn _ t'w (27rr ^ 



Qk'iJo ^x{l-x)x'{l-x') J {2txY 

(C.2) 
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with (see eq. ( ^1251) ) 

X.(^) = k-X,iq) = 75r. + ^2(1 - x) ) ' ^^''^ 

Here the variables x and x' have been introduced according to g_ = k_x 



and = k^x' . The notation is as in eq. (|3.50|) , i.e; = {si,ti,Ci). It 
is clear that for k± = the solution of this equation is of the form ( |3.5CI| ). 
By inserting this form back into ( |C.2| ) we get the eigenvalue equation for ek 
which coincides with the equation 1 + 2G n(/c) = with Il{k) given by eq. 

(PD. 



C.2 BS equation 

The BS wave function for the dynamical components of the fermion field has 
been defined in eq. (|3.47|) and satisfies the following BS equation: 



(C.4) 

where the propagator 5*++ is given in eq. ( |3.23|) , and the interaction kernel 



.25| ) has been expressed in terms of the quantity X of ( |C.3| ), where g_ = k^x 



and q'_ = k-x'. We integrate eq. (|C.4D over g+ to derive the equation for the 
wave function $^ (q) defined in eq. (^A 



^ ql +M' (fc - q)i + M- ^ ^(5). 
^ 2A;_x 2kjl-x^ 



^ -t- JKi " \ ^_ 



(27r 
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The corresponding wave function in the spin representation (see the discus- 
sion below eq. (|CTD ) 



k_J2x{l — X) 

(C.6) 

then satisfies the equation ( |C.2| ). 

D Bound state matrix elements using the LC 
BS wave functions 

If we use the BS wave function $fc in the diagram for pion decay, referring to 
the + component of the axial vector current, the resulting expression differs 
from eq. ( 2.1(j| ) only by the replacement S —>■ S for the two quark lines. 



(Note that the + component of the axial vector current can be expressed as 
^7^75'''^ •) Due to relation ( p.20| ) it is then easy to see that this replacement 



leads to no changes. If we use the wave function the trace factor involves 
two extra factors A_|_, but due to the relation 7*^7+ = \/2A^ these factors can 
be left out, leading again to the same result. 

If we use the BS wave function $fc together with the propagators S in 
the diagram fig. ^ for the pion structure function, the resulting expression 
differs from eq. (|2.13| ) only by the replacement S —>■ S for all three quark 
lines. From ( |2.2CI|) and the familiar relation S{q)'~f^'~f^ S{q) = —dS{q)/dq+ 



it follows that also 5'(g)7°7+S'(g) = —dS{q) /dq+, and performing a partial 
integration in g+ we end up with qv{x) = —g'^dll{k,x)/dk'^. Due to relation 
( |2.21D this is the same as ( [^.14| ). If we calculate the diagram fig. ^ by using 
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the wave function $fc and the propagators ( |3.23| ), the trace factor is 



Tr 



A+S{qh'^,S{q-k)A^Sl]{q-k)A+S{q-k)j'j,S{q)A+j'^^ . (D.l) 



Using the identity (see eq. ( |3.26| ) and ( p.29| ) 



S{q-k) (A+5++(g - k)A+S{q - k)) = S{q-k)A+ (1 + Z{q - k)) = S{q-k) 

(D.2) 

and 7°7^ = -\/2A+, the resuh is the same as obtained by calculating the 
diagram fig. ^ by using the propagators S. As was shown above, this gives 
the result ( p.l4| ) of the ET theory. 



^^The quantity S_^_\_ in eqs. (D.l), 
for the projection operator A+. 



is understood as the inverse of eq. (B.23) except 
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Figure captions 

1. External momenta for the qq t-matrix. 

2. Diagram for the quark LC momentum distribution in the pion. The 
double line represents the pion, the shaded area the pionic vertex func- 
tion, and cross represents the external operator 7"*". 

3. Vertex correction to the 'bare' vertex (the cross in fig. 

4. LC momentum distributions of the valence quark in the pion, using the 
LB regularization scheme. The dashed line shows the NJL model result 
(eq. (|2.26|) with the condition ( p. 36 )) with the parameters M = O.SGeV 



and A = lAlGeV . The solid line shows the result obtained by the 
QCD evolution up to next-to-leading order from Qq = O.lSGey^ to 
Q2 = AGeV^, using Aqcd = 0.25^61^ and Nj = 3. The dotted line 
shows the parametrization at = AGeV"^ obtained in ref. |]2^ by 
analyzing the experimental data. 

5. LC momentum distributions of the valence quark in the pion, using 
the TR regularization scheme. The dashed line is calculated by using 
the cut-off Atr = OASGeV for the transverse momentum integration 
in eq. ( p.26|) . For further explanation, see the caption to fig. ^. 
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